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Abstract. With the QGP opacity computed perturbatively and with the global 
entropy constraints imposed by the observed dN c h/dy 1000, radiative energy loss 
alone cannot account for the observed suppression of single non-photonic electrons. 
Collisional energy loss is comparable in magnitude to radiative loss for both light 
and heavy jets. Two aspects that significantly affect the collisional energy loss are 
examined: the role of fluctuations, and the effect of introducing a running QCD 
coupling as opposed to the fixed a s — 0.3 used previously. 

1. Introduction 

Non-photonic single electron data [TJ [2] , which present an indirect probe of heavy quark 
energy loss, have significantly challenged the underlying assumptions of jet tomography 
theory. A much larger suppression of electrons than predicted [3] was observed in the 
p T 4 — 8 GeV region. "These data falsify the assumption that heavy quark quenching 
is dominated by [pQCD based] radiative energy loss when the bulk [weakly coupled] 
QCD matter parton density is constrained by the observed dN/dy ~ 1000 rapidity 
density of produced hadrons." [4] 

WHDG [1] revisited the assumption that pQCD collisional energy loss is negligible 
compared to radiative energy loss [HIE]. As argued there, and references therein, "the 
elastic component of the energy loss cannot be neglected when considering pQCD jet 
quenching." As shown in WHDG and elsewhere [7j, the computationally expensive 
integrations over the geometry of the QGP cannot be reduced to a simple 'average 
length' prescription. Indeed, this computation time is essential to produce radiative + 
collisional energy loss calculations consistent with the pion data. 

There are large theoretical uncertainties in the WHDG results [8] . Very significant 
to the electron prediction is the uncertainty in the charm and bottom cross-sections. 
There are also theoretical uncertainties in the energy loss mechanisms. Here, two aspects 
of the collisional energy loss will be examined with the aim of improving the energy loss 
model. 

2. Collisional energy loss fluctuations 

Similar to radiative energy loss, the fluctuations of collisional energy loss around the 
mean affect the quenching of the quark spectra. Collisional fluctuations are often 
modelled in a Fokker-Planck formalism, characterized by two numbers or functions: 
drag and diffusion. WHDG implemented an approximation to this scheme applicable 
for small energy loss by giving the collisional loss a gaussian width around the mean, 
with a 2 = 2T(e), where (e) is the mean energy loss given by a leading log calculation. 
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The drag-diffusion method is essentially a continuum approximation to a discrete 
process. A high energy jet traversing the QGP will undergo only a small number of 
collisions. In the Gyulassy-Wang model, the expected mean free path of a quark is 
~ 2fm, so there is a very significant surface region in which the fluctuations will differ 
greatly from those given by the continuum approximation. It is therefore necessary to 
look at the fluctuations per collision and in the number of collisions. A simple model 
to investigate this is to model the medium as initially static objects which will then 
recoil upon collision, model the interaction between jet and medium using the full HTL 
medium modified propagator. This gives the probability of longitudinal momentum loss: 
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This single collision distribution is then Poisson convoluted to give the distribution 
for a finite number of expected collisions: 

e -(N coll {E)) . 

P„+i(e, E)— — — — / dx 1 ... dx n p(x 1 , E) . . . p(x n , E)p(e-x 1 -. . -x n , E)(2) 
[n+ly. J 

The mass of the medium particle is tuned to give an average energy loss similar to 
that of the BT and TG leading log calculations (m ~ 0.2GeV - although here we are 
not interested in the average energy loss per se). In Fig. [TJ the probabiliy of fractional 
energy loss in one collision is shown, similar to a 1/t 2 Bjorken collisional style model, 
with screening at small t- values (shown in the right pane of Fig. [[]). 
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Figure 1. The distribution of fractional energy loss e (left) and Mandalstam variable 
t (right - scaled by t 2 ) in a collision using this model. 

Figure [2] illustrates the distributions in energy loss for a finite number of collisions 

for bottom and light quark jets. The results for charm quarks are qualitatively similar to 

those for light quarks. For a large number of collisions (eg average number of collisions 

(n) = 10, L~ 20fm), the distributions are roughly symmetric and somewhat similar to 

the simple WHDG gaussian. This is expected from the central limit theorem. The Raa 

values extracted from these distributions are similar, with (n) = 10 and the gaussian 

approximation only differing by ~ 0.01. Surprisingly, a similar result for the Raa 

values is found for (n) = 2 collisions for bottom quarks. The large change arrives 

for light quarks. For both (n) = 2,5 collisions, the gaussian approximation gives a 

very different distribution for the fluctuations and a very different Raa value. The 
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Figure 2. The distribution in fractional energy loss for a 20 GeV jet, for a bottom 
quark jet and a light quark jet, for different numbers of collisions. Shown are lines for 
the gaussian approximation, exactly n collisions (labelled eg 'n=10') and on average 
n collisions (labelled eg '(n) = 10'). Inset is the Raa for these distributions and the 
Raa evaluated for a delta function distribution at the mean loss (+ point). 

gaussian approximation overpredicts the Raa suppression by 0.1, which is around a 
30% effect for (n) = 2 collisions. This cannot be neglected. A full treatment of the finite 
number of collisions will reduce the quenching due to elastic energy loss compared to the 
treatment in WHDG. This conclusion is also applicable to other uses of Fokker-Planck 
/ Langevin formalisms that use a continuum description of the collisional process. The 
Raa predictions for bottom quarks are likely only marginally affected, those for light 
quarks most affected. 
3. Running QCD coupling 

In [8], the change of the fixed QCD coupling a s from 0.3 to 0.4 was seen to significantly 
change the Raa precitions from the WHDG model. There has been much recent work 
on the effect of a running coupling on the collisional energy loss 0, [10] (ie a s = a s {Q 2 )). 
Here, we revisit the collisional energy loss in a similar manner to [9] , looking at a simple 
Bjorken-style estimate [llj. Bjorken's estimate for the collisional energy loss is: 

^-(§f*(^)<*M?) 

In [H], the running coupling version for very high jet energies is given as: 

although this neglects the finite energy kinematic bound on the jet. Adding in this 
bound to this calculation gives 

dE n n / 1 \ . „o „ „. . „o n.„n, / ATE 
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which is similar in structure to the original fixed coupling estimate. A numerical 
comparison of equations l3"]T4"fo1 is shown in Fig. [21 For reasonable temperatures, 
T ~ 0.2 — 0.3GeV, all results are of a similar order of magnitude. For reasonable 
energies, no qualitatively new behavior is seen (although, as found in [9], the E — » oo 
behavior is new, but this affects much higher energies than those of interest at RHIC 
or even LHC). When the kinematic bounds are taken into account, the result for the 
average energy loss including running coupling is often larger than the fixed a s = 0.3 
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result used in [3]. However, the numerical result is very sensitive to the input parameters 
used, illustrated in the middle and right panes of Fig. [3] for changing the prescription 
for /i from that used in |3] to that from [9]. 
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Figure 3. Left: the average energy loss for a light quark jet evaluated for (1) 
Peshier's running a s — a s (Q 2 ) but infinite energy jet approximation (2) Finite energy 
running a s , (3) fixed a s independent of Q 2 but evaluated at Q 2 = (2ttT) 2 and (4) fixed 
a s = 0.3. The functional form of a s = a s (Q 2 ) is taken for vacuum at 1-loop as in [9]. 
Middle and right: the ratio between the different versions for two different evaluations 
of the Debye screening scale fi, (1) 'self-consistent' /z[9] and (2) fi evaluated at a fixed 
temperature. 'Bjorken (2)' is with all parameters evaluated with a s = 0.3. 

4. Conclusion 

It has been argued previously "that radiative and elastic average energy losses for heavy 
quarks were in fact comparable over a very wide kinematic range" [3] , and even "E > 
10 GeV light and charm quark jets have elastic energy losses smaller but of the same 
order of magnitude as the inelastic losses" |4J. Hence, collisional energy loss cannot be 
neglected when considering jet quenching of high energy jets in the QGP at either RHIC 
or LHC. A simple model combining collisional and radiative energy losses significantly 
reduces the discrepancy between the predictions and data. 

Two possible improvements to the WHDG model have been examined here. The 
inclusion of a finite number of collisions is seen to reduce the effect of the collisional 
energy loss on the quenching of gluons, light and charm quarks, but not to significantly 
affect the bottom quark Raa- Opposite to this effect, including a running QCD coupling 
increases the energy loss by up to a factor of 1.5. The combination of these two affects, 
along with other large uncertainties in the prediction for electron Raa such as the ratio 
of charm to bottom total cross-sections, hints at the possibility that both the pion and 
electron i^^^s may both be within range of purely perturbative calculations. 
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